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Probing spectral properties of radio-quiet quasars searched 
for optical microvariability 
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ABSTRACT 

We obtained SDSS spectra for a set of 37 radio-quiet quasars (RQQSOs) that had 
been previously examined for rapid small scale optical variations, or microvariability. 
Their H/3 and Mg II emission lines were carefully fit to determine line widths (FWHM) 
as well as equivalent widths (EW) due to the broad emission line components. The 
line widths were used to estimate black hole masses and Eddington ratios, t. Both 
EW and FWHM are anticorrelated with £. The EW distributions provide no evidence 
for the hypothesis that a weak jet component in the RQQSOs is responsible for their 
microvariability. 
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1 INTRODUCTION 

Over the past 15 years there have been rather extensive 
examinations of a significant sample of radio-quiet QSOs 
(RQQSOs) and Seyfert galaxies for small brightness changes 
(typically 0.02 mag) over short times (a few hours) (e.g., 
Stalin et al. 2004b; Carini et al. 2007; Ramirez et al. 2009). 
This phenomenon of microvariability, or intranight opti- 
cal variability (INOV) was first confirmed for blazars (e.g., 
Miller, Carini & Goodrich 1989; Carini 1990) for which mi- 
crovariability almost certainly arises from the relativistic jet 
(e.g., Marsher, Gear & Travis 1992), at least when the source 
is in an active state; however, in low states it is possible 
that rapid fluctuations are due to processes originating in or 
just above the accretion disc (for a review, see Wiita 2006) . 
Since RQQSOs lack significant jets, the microvariability in 
radio-quiet objects may arise from processes on the accre- 
tion disc itself, and thus could possibly be used to probe the 
discs (e.g., Gopal-Krishna, Wiita & Altieri 1993; Stalin et 
al. 2004a; Kelly, Bechtold & Siemiginowska 2009). 

Recently Carini et al. (2007) reported new observations 
for several sources and also compiled a sample of 117 objects 
from the literature which have been searched for microvari- 
ability (their Table 3). Of these, 47 are classified as Seyfert 
galaxies, 6 as broad absorption line (BAL) QSOs, and 64 as 
QSOs. In addition, in order to learn which classes might be 
the most fruitful for making further searchs for microvari- 
ability they also classified their sample in terms of following: 
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redshift distribution; radio loudness, through R, the ratio of 
the radio [5 GHz] flux to the optical [4400 A] flux); optical 
magnitude, my; luminosity; and observing strategy. In their 
entire sample 21.4% of the objects were found to exhibit mi- 
crovariability, but among objects classified as Seyfert galax- 
ies, BAL QSOs and QSOs, microvariability was found in 
17%, 50% and 23.2%, respectively (Carini et al. 2007). The 
observed high fraction of microvariations in BAL QSOs (al- 
though the sample is quite small) suggests that while plan- 
ning microvariability studies to investigate the physical pro- 
cesses in or near the accretion disc, it might be worthwhile 
to invest more observing time on the BALQSO class. 

With improvements in observation quality as well as 
in sample size, constraints on the models capable of pro- 
ducing the intranight variability have also improved over 
last decade. Recently, Czerny et al. (2008) have used non- 
simultaneous optical and X-ray data of 10 RQQSOs with 
confirmed INOV to compare observational constraints on 
the variability properties with the predictions of theoreti- 
cal models such as: (i) irradiation of an accretion disc by 
a variable X-ray flux (e.g., Rokaki, Collin- Souffrin & Mag- 
nan 1993; Gaskell 2006); (ii) an accretion disc instability 
(e.g., Mangalam & Wiita 1993); (iii) the presence of a weak 
Doppler boosted jet, or "blazar component" (e.g., Gopal- 
Krishna et al. 2003). Their investigation suggests that a 
blazar component model yields the highest probability of 
detecting INOV. 

In this blazar component scenario, spectral properties 
of the sources can play a crucial role in constraining the 
models further. For instance, if blazar components are dom- 
inating the variability of RQQSOs, then, due to the increase 



2 Chand, Wiita & Gupta 



in the continuum level, one would expect emission lines to 
be diluted. Therefore smaller equivalent widths (EWs) of 
prominent emission lines such as H/3 and Mg II should be 
detected in sources that showed microvariability when com- 
pared to their average values in the whole sample. If we take 
the extreme case of BL Lacertae objects, which often lack 
observable emission lines and are usually defined as objects 
that have no emission line with an EW ^5 A (e.g., Stickel 
et al. 1991; cf., Marcha, et al. 1996), this dilution by the jet 
component is understood to be severe. So it becomes very 
important to test whether microvariability of RQQSOs has 
any correlation with spectral parameters such as EW and 
full width at half maximum (FWHM) of prominent emission 
lines. The H/3 and Mg II emission lines are very promising for 
such investigations, as these lines have also been found to be 
very useful in estimating other key parameters of AGN cen- 
tral engines such as black hole (BH) mass (e.g., McGill et al. 
2008) and Eddington ratio (e.g., Dong et al. 2009a, 2009b). 
The average H/3 EW of a large sample of quasars is found 
to be around 62. 4A (Foster et al. 2001), so any correlation 
of H/? EW or FWHM would not only give insight about 
the nature of variability, such as the presence or absence of 
blazar components, but also would be very useful for mak- 
ing a promising sample for future microvariability studies. 
Similarly the measurements of BH masses for RQQSOs that 
show microvariability may well be another important con- 
straint on the models trying to understand the nature of 
their optical microvariability. 

Here we have worked toward these goals by exploiting 
the optical spectra available from Sloan Digital Sky Survey 
(SDSS) Data Release 7 (DR7; Abazajian et al. 2009) with 
careful spectral modeling of the H/3 and Mg II emission line 
regions. First we aim to investigate any effect of these key 
spectral parameters (e.g., EW and FWHM) on microvari- 
ability of RQQSOs. Second we estimate other relevant AGN 
parameters such as the black hole mass, Mth, and the Ed- 
dington ratio in the context of our RQQSOs sample, which 
has been extensively searched for microvariability. 

The paper is organized as follows. Section 2 describes 
the data sample and selection criteria while Section 3 gives 
details of our spectral fitting procedure. In Section 4 we fo- 
cus on BH mass measurements and in Section 5 we give 
estimates of Eddington ratios and of BH growth times. Sec- 
tion 6 gives a discussion and conclusions. Throughout, we 
have used a cosmology with Hq=70 kms -1 Mpc -1 , fiM=0.3 
and fU=0.7. 



2 DATA SAMPLE AND SELECTION 
CRITERIA 

Our sample is chosen from the compilation of 117 radio- 
quiet AGNs in Carini et al. (2007). These sources have all 
been searched, often extensively, for optical microvariabil- 
ity. When there have been several attempts at monitoring 
the same RQQSO or Seyfert galaxy by different groups of 
authors around the globe they often have used different pop- 
ular names or even different IAU names for the same source. 
Therefore we first looked in all the papers cited in the com- 
piled data (Carini et al. 2007, their Table 3) and obtained the 
right ascensions (02000. o) and declinations(520oo.o) of all the 
sources monitored by different groups involved in searches 



for their optical microvariability. Then we searched these 
sources in the recent catalog of AGN and Quasars by Veron- 
Cetty & Veron (2006) and made the list of all sources mon- 
itored until 2005 by giving their popular name from Veron- 
Cetty & Veron (2006). 

We searched for consistently produced optical spectra 
of these sources by looking for these targets in SDSS DR7 
(Abazajian et al. 2009) and found that out of these 117 
sources (Carini et al. 2007) SDSS DR7 spectra were avail- 
able for the 37 sources listed in Table [1] The first four 
columns give identification information and redshift, the 
next two provide apparent and absolute V magnitudes and 
the last two indicate whether microvariability was detected 
and whether the radio quiet AGN has been classified as a 
QSO, BALQSO, Seyfert (Sy) galaxy, or Narrow Line Seyfert 
1 (NLS1) galaxy. From among the 37 sources in our sample, 8 
source spectra cover both Mg n and H/3 lines, 11 spectra cov- 
ered only H/3 lines while 18 spectra covered only Mg 11 dou- 
blet lines. This set of RQQSOs becomes an unique sample 
to look for any dependence on microvariability of their op- 
tical spectral properties as its members had been already 
targeted to search for optical microvariability in the past. 



3 LINE ANALYSIS THROUGH 

SIMULTANEOUS SPECTRAL FITTING 

The spectra were first corrected for Galactic extinction using 
the extinction map of Schlegel et al. (1998) and the redden- 
ing curve of Fitzpatrick (1999). Then they were transformed 
into the rest frame using the redshift given in the header of 
the SDSS spectra. 

In order to have accurate measurements of spectral 
properties such as the EWs and FWHMs of emission lines, 
one has to carefully take into account the contribution from 
continuum as well as from Fe 11 emission multiplets. The 
common practice was to first fit and subtract the AGN con- 
tinuum and then the Fe 11 emission (e.g., Boroson & Green 
1992; Marziani et al. 2003) and finally the other emission 
lines. The other possibility is to fit simultaneously the con- 
tinuum and the Fe n emission multiplets and other emission 
lines (e.g., Dong et al. 2008, 2009a). Unfortunately, for the 
optical spectra of most of RQQSOs in our sample, fitting 
the continuum and the Fe 11 emission is very complicated 
because: (i) there are essentially no emission-line free re- 
gions where the continuum can be determined (Vanden Berk 
et al. 2001); (ii) the prominence of Fe n features and their 
blending with the H/3 and Mg n lines; (in) the H/3 line is 
highly blended with the O inAA4959,5007 lines. As a result, 
we have opted to carry out simultaneous fit£| of continuum, 
Fe n emission, H/3 and Mg 11 and all other metal emission 
lines present in the spectra. 

For this purpose we follow the following procedure. 



1 To carry out the simultaneous fit we have used the 
MPFIT package for nonlinear fitting, written in In- 
teractive Data Language routines. MPFIT is kindly 
provided by Craig B. Markwardt and is available at 
|http://cow.physics.wisc.edu/~craigm/idl/| 
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Table 1. Our sample of radio-quiet QSOs and Seyfert galaxies from the compilation of Carini et al. (2007) 
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3.1 H/3 region fit 

For fitting the optical region comprising the H/3 line, we con- 
sider the rest wavelength range between 4435 and 5535 A. 
The continuum is modeled by a single power law, as since we 
limit ourselves to below 5600A there is no need of a double 
or broken power law fit (e.g, Dong et al. 2008). As the profile 
of the H/3 line is rather complex, it is fitted with multiple (as 
many as four) Gaussians. As an initial guess we consider two 
narrow, one broad and one very broad components. To re- 
duce the arbitrariness of the component fits and also to make 
the decomposition more physical, we have constrained the 
redshift and width of the two narrow components of H/3 to 
be the same as those of the O mAA4959,5007 lines. The line 
profile of O IIIA5007 (and hence of O IIIA4959 ) are modeled 
as double Gaussians, with one stronger narrow component 
with width less than 2000 km s _1 and one weaker broader 
component with width less than 4000 km s _1 . However, if 
for some source spectra the second component is not statis- 



tically required, the procedure we use automatically drops 
it during the fit. So, basically we have only four free pa- 
rameters for the width and redshift of O mAA4959,5007 and 
narrow H/3 components (two for redshift and two for width). 

The optical Fe 11 emission is modeled with two separate 
sets of templates in analytical forms similar to those used by 
Dong et al. (2008). One template is for the broad-line sys- 
tem and the other for the narrow-line system, and we assume 
the form (7(A) = Cf,C(,(A) + c n C„(A), where Ct,(A) represents 
the broad Fe n lines and C„(A) the narrow Fe 11 lines, with 
the relative intensities fixed at those of IZW1, as given in 
Tables A.l and A. 2 of Veron-Cetty, Joly & Veron (2004). 
The redshift of the broad Fe 11 lines and the H/3 component 
is fitted as a free parameter while their widths are kept the 
same, with a constraint that they should be larger than 1000 
km s _1 ; this width is used for estimating the BH masses. 
Similarly for the narrow Fe 11 line the redshift is fitted as 
a free parameter but the width is kept the same as that 
of the stronger narrow O IIIA5007 component. The fourth 
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Figure 1. The best fit to the H/3 emission line of the SDSS spectra of the QSO CSO 448. Upper panel: complete spectrum fit (thick 
dashed/red) and components of the fit: power law continuum (thin dashed/black line), broad Fe II (dot-dashed/blue), narrow Fe II (dot- 
ted/magenta) lines. Lower panel: continuum, Fe II and metal line subtracted spectrum (solid/black) with the best fit total H/3 profile 
(thick dashed/red) and H/3 components (dotted/blue) lines. Note that the entire fit is performed simultaneously (not first continuum 
subtraction then H/3 fit) but these aspects are shown separately for the sake of clarity. 



(very broad) H/3 component, if required for the fit, is sub- 
ject only to the constraint that its width should be more 
than 1000 km s _1 . The emission lines other than Fe II and 
O mAA4959,5007 (see Table 2 in Vanden Berk et al. 2001), 
are modeled with single Gaussians. The final fit is achieved 
by simultaneously varying all the free parameters to mini- 
mize the x 2 value until the reduced Xr is ~ 1- Samples of 
our spectral fitting in the optical region are given for H/3 in 
Fig. [T] and Fig. [2] The values for FWHM and EW (both the 
broad component, EWb, and total, EW„n) for the H/3 lines 
are given in Table [2] 



3.2 Mg II doublet region fit 

For fitting the Mg II doublet region, we consider a rest wave- 
length range between 2200 and 3200 A. The continuum in 
this region is modeled by a power law, i.e., aiA _c \ For fit- 
ting the Mg n doublet we have tried both Gaussians and 
truncated 5-parameter Gauss-Hermite series (Salviander et 
al. 2007; Dong et al. 2009a). In Gaussian profile fits, our ini- 
tial guess consists of two Gaussian components for each line 
of the Mg II doublet; however, whenever the second com- 
ponent is not statistically required the procedure automat- 



ically drops it during the fit. Only two cases where single 
components fits were found to be sufficient were detected 
(UM 479 and US 1443). To test for any overfitting caused 
by assuming two components, we also forced our procedure 
to fit only single components; however, in doing so for all 
other sources, a very good fit is found neither for the wings 
of the lines nor for their central narrow cores. The trun- 
cated 5-parameter Gauss-Hermite series was also tried, but 
the fits in this case also were not as good as those from the 
two component Gaussian profile model, so our final analyses 
employ that fitting approach. 

In our two component Gaussian profile fits, the red- 
shift and width of each component (narrow/broad) of 
Mg IIA2796 were tied to the respective components of the 
Mg IIA2803 line. In addition, we have constrained the width 
of narrow component to be smaller than 1000 km s _1 and 
the width of the broader Mg II component to be same as 
width of UV Fe II emission line in the region (UV Fe n ). 
We used an UV Fe II template generated by Tsuzuki et al. 
(2006), basically from the measurements of IZW1, which 
also employ calculations with the CLOUDY photoioniza- 
tion code (Ferland et al. 1998). This template is scaled and 
convolved to the FWHM value equivalent to the broad com- 
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Figure 2. Same as Fig. [T] for CSO 879. 



Table 2 . Results of the H/3 line analyses 



QSO Name 
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L(5100A) 


FWHM(km s" 1 ) 




Lbol 


EW B (A) 


EW„n(A) 




10 44 crg s- 1 


L edd 
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0.07 
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10461.15 


9.42 


0.02 


48.81 


75.77 + 0.49 


PG 0931+437 





458 


31.42 


9537.59 


9.62 


0.05 


11.47 


68.42 + 0.52 


PG 1049-005 





359 


24.17 


5193.56 


9.03 


0.15 


22.16 


76.06 + 0.40 


PG 1259+593 





474 


42.13 


2922.88 


8.65 


0.63 


4.33 


41.06 + 0.50 


PG 1307+085 





154 


5.22 


4080.08 


8.49 


0.11 


62.77 


99.83 + 0.54 


PG 1309+355 





183 


10.98 


4326.57 


8.70 


0.15 


20.86 


35.27 + 0.46 


PG 1444+407 





268 


10.15 


3528.12 


8.51 


0.21 


8.26 


66.90 + 0.50 


Q 1252+0200 





344 


16.62 


23694.59 


10.27 


0.01 


26.78 


90.09 ± 0.54 


TON 52 





434 


16.68 


2724.97 


8.39 


0.46 


11.90 


38.78 + 0.69 


US 1867 





514 


25.11 


2837.29 


8.52 


0.52 


33.89 


68.80 ± 0.46 


US 3150 





469 


14.44 


3402.68 


8.55 


0.27 


14.78 


68.28 + 0.67 


US 3472 





532 


25.48 


4587.36 


8.94 


0.20 


14.30 


69.93 ± 0.64 


US 995 





227 


3.90 


4509.77 


8.51 


0.08 


46.58 


51.95 + 0.58 
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Figure 3. Mg II emission line fits for two RQQSOs, CSO 21 and CSO 18. The upper plot in each panel shows the SDSS spectra and 
continuum as a dashed line with the best fit UV Fe II template as a solid (magenta) line. The lower plot in each panel show the continuum 
and Fe II subtracted spectrum (solid) the best fit Mg II profile (thick dashed/red) and Mg llAA2796,2803 components by (dotted/blue) 
lines. Note that the entire fits are performed simultaneously (not first continuum subtraction, then Mg II fit) but are shown separately 
for the sake of clarity. 



ponents of Mg II by taking into account the FWHM of the 
I ZW 1 template. The best fit value of the broad component 
of Mg n obtained in this way is finally used in our calcu- 
lation of BH mass. The emission lines other than Fe n lines 
identified from the composite SDSS QSO spectrum (see Ta- 
ble 2 in Vanden Berk et al. 2001), are modeled with single 
Gaussians. The final fit is achieved by varying all the free 
parameters simultaneously by minimizing the \ 2 value, un- 
til the reduced Xr is ~ 1- Demonstrations of our spectral 
fitting in the UV region are given in Fig. [3] and values for 
FWMHs and EWs are given in Table [3] 



3.3 Optical microvariability and spectral 
properties 

Among the 19 sources for which we have spectral coverage 
of the H/3 line, optical microvariability is shown by 7 of them 
while the other 12 have not been seen to show this prop- 
erty (Carini et al. 2007). Similarly, among the 26 sources 
with spectral coverage of Mg n doublet, only 5 show optical 
microvariability but the other 21 do not. 

We have used our best fit values of FWHM and EW 
to search for any correlation between them and optical mi- 



crovariability. The median value of total rest frame EW of 
the H/3 line (i.e., the sum of the EWs of all fitted com- 
ponents) for 12 sources without optical microvariability is 
found to be ~ 69 A. Its value for the 7 sources having opti- 
cal microvariability is ~ 79 A. This seems to be contrary to 
the expectation based on a blazar component model of mi- 
crovariability for RQQSOs, which would predict smaller EW 
values for variable sources, due to dilution of emission line 
strength by jet components (Czerny et al. 2008). However, it 
should be noted that observed H/3 profiles certainly require 
multiple components to fit them. Therefore, if optical mi- 
crovariability is a phenomenon related to the central AGN 
engine, then it is presumably better to do such a comparison 
only using the broad component of the H/3 line (Section l3.1[l . 
This is because the clouds responsible for this broad com- 
ponent are clearly in the sphere of influence of the massive 
BH, while the other components may not be. The median 
value of rest frame EW of this component for objects show- 
ing any optical microvariability is 26.8 A and for objects 
without optical microvariability it is 20.8 A, which is again 
contrary to the expectation based on a blazar component 
model of optical microvariability. The corresponding values 
of median FWHM of these broad components for sources 



Spectra of microvariable radio-quiet QSOs 7 




HjSEW [A] Mgll EW [A] 



Figure 4. Histograms of FWHM and rest frame equivalent width (EW) based on best fits of H/3 and Mg II lines. The shaded regions 
correspond to sources with confirmed optical microvariability while the non-shaded regions are for those sources for which optical 
microvariability is not detected. 



with and without microvariability are somewhat different, 
at 4509 km s _1 and 4326 km s _1 , respectively. 

Similar to what we did with the H/3 line, we have also 
used our Mg II spectral fits to look for any correlation with 
microvariability. The median value of rest frame EW of 
broad Mg II doublet components (with width tied to the 
UV Fe II template width) is found to be 9.15 A for the five 
sources with optical microvariability and 9.24 A for the 21 
sources without any confirmed microvariability. These val- 
ues are essentially identical, and so again do not support 
the simple blazar component model for microvariability in 
RQQSOs. The corresponding values of median FWHM for 
sources with and without microvariability are 2432 km s _1 
and 2447 km s _1 , respectively, and so are indistinguishable. 

In Fig. [4] we show the histograms of our FWHM and 
rest frame EW values based on our best fits for H/3 and 
Mg II lines. The shaded and non-shaded regions correspond 
respectively to sources with and without confirmed optical 
microvariability. As the figure shows, the distributions of 
sources with and without microvariability are on the whole 
quite similar, though there seems to be a hint that the 
median EW width of the H/3 line might be larger for 
sources with optical microvariability compared to 
those without optical microvariability. To quantify 
this posibility we performed Kolmogorov-Smirnov 



tests on all the distributions shown in Fig. [31 The 
probabilities that the sources with and without 
confirmed optical microvariability are drawn from 
similar distributions are: 0.613 for H/3(EW); 0.664 
for Mg II(EW); 0.792 for H/3(FWHM) and 0.995 for 
Mg II(FWHM). As these null probabilities are all high 
our data samples do not establish or support any re- 
lation of FWHM or EW with the optical microvari- 
ability properties. The lack of a difference in the EW 
values does not support simple models in which the 
microvariability arises from jets. 



4 BLACK HOLE MASS MEASUREMENTS 

The strong Mbh-<? correlation between supermassive black 
hole (SMBH) mass and the velocity dispersion of the sur- 
rounding spheroidal stellar system (e.g., Tremaine et al. 
2002) have shown that black hole (BH) assembly is inti- 
mately related to the bulge properties of its host galaxy. 
This paradigm of "BH-galaxy co-evolution" boosted the ef- 
forts to measure BH masses over a wide range (e.g., Rich- 
stone et al. 1998; Nelson 2000; Bower et al. 2000). However, 
one cannot obtain BH mass estimates from either stellar 
or gas dynamics, based upon traditional methods, beyond 
~100 Mpc, due to the difficulty in resolving the BH sphere 
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Table 3. Results of the Mg II line analyses 



QSOName z ^"""^i FWHM(kms~ 1 ) Log(^£ Log(^£ T± EW B (A) EW aII (A) 



CSO 18 




1.300 


139.87 


2728.10 ±0.49 


8.71 


9.25 


0.32 


20.28 


22.09 ±0.37 


CSO 21 




1.190 


151.63 


3244.49 ± 0.82 


8.88 


9.41 


0.24 


22.80 


24.65 ±0.33 


CSO 233 




2.030 


220.37 


1670.96 ± 0.95 


8.40 


8.92 


1.08 


17.60 


17.60 ±0.84 


CSO 879 




0.549 


44.02 


1548.88 ±0.50 


7.90(8.85) 


8.50 


0.56 


10.47 


14.96 ±0.33 


PG 0935-1 


-416 


1.966 


711.91 


2518.34 ±1.19 


9.08 


9.53 


0.86 


4.68 


5.00 ±0.34 


PG 0946-1 


-301 


1.220 


355.37 


2216.44 ±0.37 


8.78 


9.27 


0.78 


13.62 


14.80 ±0.27 


PG 1206n 


-459 


1.155 


607.67 


2441.60 ±0.76 


9.01 


9.47 


0.84 


12.30 


13.46 ±0.21 


PG 1248- 


-401 


1.032 


254.23 


2337.98 ±0.28 


8.73 


9.24 


0.59 


16.82 


18.84 ±0.30 


PG 1254- 


-047 


1.018 


231.12 


3603.36 ±0.57 


9.08 


9.60 


0.24 


15.65 


18.63 ±0.28 


PG 1259n 


-593 


0.472 


76.12 


2417.71 ±0.65 


8.44(8.65) 


9.01 


0.30 


8.77 


12.71 ±0.28 


PG 1338n 


-416 


1.204 


249.35 


3062.60 ± 0.85 


8.96 


9.47 


0.34 


13.22 


16.08 ±0.29 


PG 1630-1 


-377 


1.478 


612.30 


2793.42 ± 1.91 


9.13 


9.59 


0.65 


12.12 


13.45 ± 0.26 


Q 1628.5+3808 


1.461 


193.09 


2167.82 ±0.58 


8.59 


9.12 


0.60 


16.54 


16.54 ±0.37 


Q J07514 


2919 


0.912 


347.62 


2432.96 ± 0.37 


8.85 


9.34 


0.64 


14.79 


15.98 ±0.25 


UM 497 




2.022 


342.38 


1901.77 ±0.56 


8.63 


9.13 


1.04 


22.83 


22.83 ±0.46 


US 1420 




1.473 


199.32 


1708.34 ± 0.45 


8.40 


8.92 


0.98 


7.50 


10.76 ±0.33 


US 1443 




1.564 


302.09 


2589.39 ± 0.61 


8.87 


9.37 


0.53 


14.10 


14.10 ±0.30 


US 1498 




1.406 


162.70 


1844.98 ±0.29 


8.41 


8.94 


0.76 


11.73 


12.87 ±0.34 


US 1867 




0.513 


39.58 


1745.21 ±0.39 


7.98(8.52) 


8.58 


0.42 


13.12 


16.87 ±0.24 


US 3150 




0.467 


22.81 


1467.23 ± 1.09 


7.68(8.55) 


8.31 


0.45 


8.84 


13.46 ±0.39 


US 3472 




0.532 


45.26 


2234.99 ± 0.39 


8.23(8.94) 


8.83 


0.27 


15.13 


18.94 ±0.32 


TON 52 




0.434 


28.10 


1739.86 ±0.78 


7.88(8.39) 


8.51 


0.36 


17.31 


17.31 ±0.49 


CSO 174 




0.653 


36.73 


3564.17 ±1.06 


8.58(9.24) 


9.19 


0.10 


32.41 


36.94 ±0.36 


TON 34 




1.925 


1761.85 


3006.16 ±0.62 


9.47 


9.88 


0.94 


10.77 


14.04 ±0.28 


CSO 179 




0.782 


95.95 


2691.86 ±0.28 


8.59(9.18) 


9.15 


0.28 


20.44 


25.62 ±0.33 


PG 1522- 


-101 


1.328 


765.86 


3437.96 ± 0.90 


9.37 


9.82 


0.48 


15.66 


16.13 ±0.25 



a Using the McLure & Dunlop (2004) scaling relation, i.e., Eq.[2] The values given in parentheses are from 
our H/3 measurements for the same object from Table [2] 

k Using the fixed slope of the i — L relation from Dietrich et al. (2009), i.e., Eq. [3] 



of influence. For these more distant sources the only viable 
way to probe BH masses is through their AGN activity. 

The reverberation mapping technique (e.g., Blandford 
& McKee 1982; Peterson 1993) is the most direct method 
to probe gas dynamics on spatial scales close to the BH in 
AGN, but the required long-duration spectroscopic moni- 
toring campaigns currently preclude the method from being 
applied to large numbers of objects. Because the time-lags 
involved grow linearly with the SMBH mass, only relatively 
nearby Seyfert galaxies with modest values of Mbh have 
been studied in this fashion. Therefore for any large statisti- 
cal studies one has to rely on the more indirect technique of 
virial single-epoch methods, which were pioneered by Dibai 
(1980) and shown to be consistent with reverberation map- 
ping masses (e.g., Bochkarev & Gaskell 2009). These mea- 
surements make use of the widths of emission lines such as 
H/3 , Mg n or C IV in conjunction with the bolometric lu- 
minosity (e.g., Vestergaard & Peterson 2006). This method 
relies upon the FWHM indicating the typical velocities of 
clouds bound to the BH and the assumption that the mean 
distances, r, of the clouds emitting particular lines are de- 
termined by their ionization potentials, and thus scale with 
the ionizing continuum, L, of the QSO. We have used our 
H/3 or Mg II FWHM measurements to compute the BH mass 
of RQQSOs in our sample, as such estimates will also im- 
pose important constraints in studies trying to understand 
the nature of RQQSOs (e.g., Czerny et al. 2008). 



Vestergaard & Peterson (2006) presented improved em- 
pirical relationships useful for estimating the central BH 
mass based on FWHMs of emission lines such as H/3 , 
Mg II and C IV . We used their following scaling relation 
to estimate the BH masses based on H/3 lines (Vestergaard 
et al. 2006, their Eq. 5) 



log Af BH (H/3) = log 



( FWHM(H/3) V 
^1000 km s- 1 J 

+ (6.91 ± 0.02) + log f ^5100A) \ 
v V 10 44 erg s- 1 J 



0.50±0.06 



(1) 



where La (5100 A) is the monochromatic luminosity at 
5100 A, which we have computed from the best fit power- 
law continuum, a\\~ a , in our simultaneous fit of the whole 
spectral region (see Section [23}. 

As is evident from this equation, careful measurements 
of the FWHM of emission lines are crucial for good BH mass 
estimations. So it is very important to decompose the lines 
into broad and narrow components, since the former are 
much more likely to be in the BH sphere of influence, while 
the later are contributed from the much more distant nar- 
row line region under the dynamical influence of the galaxy's 
stellar component. Therefore, as described in Section r3.ll we 
have constrained the narrow component's profile to be linked 
to the O iiiAA4959,5007 line profile while the broad compo- 
nent's profile is associated with the broad Fe II emission mul- 
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Figure 5. Comparison of BH masses estimated using H/3 and 
Mg II lines for the eight sources for which spectra cover both 
H/3 and Mg II lines. The Mg II based BH masses using Eq. [2] 
(McLure & Dunlop 2004) are shown by hollow circles, while those 
based on Eq. [3] (Dietrich et al. 2009) are shown by filled circles. 
The dotted line shows equal BH masses estimated using H/3 and 
Mg II lines. 



tiplet. We have modeled the whole Fe II emission multiplet 
by also decomposing it into broad and narrow component 
contributions. This aspect of the fitting has been largely 
ignored in almost all previous BH mass estimation studies 
until recent efforts by Dong et al. (2008, 2009a). So in the 
above Eq.[T]the FWHM we used, from among our Gaussian 
de-composition of H/3 line (Section [3}, is the one with its 
width tied with the width of the broad Fe II emission line, 
after correcting for instrumental broadening. The result of 
our analysis using H/3 lines is summarized in Table [2] 

To estimate black hole masses based on Mg n , we ini- 
tially adopted the calibrations in McLure & Dunlop (2004) , 



log Mbh (Mg II) = (0.62 ± 0.14) log 



A£ A (3000A) 
10 44 erg s" 1 



+ 21og ^F\VllM(M,^ II) 



km s~ 



+ 0.505 



(2) 



where the monochromatic luminosity at 3000 A has been 
computed from our power-law fit of the continuum and the 
FWHM used is the broad component of our double Gaussian 
decomposition model, after correcting for the instrumental 
broadening. The result of our Mg II analysis for black hole 
mass based on the McLure & Dunlop (2004) scaling relation 
is summarized in Column 5 of Table [3] 

In our sample eight QSO spectra cover both Mg II and 
H/3 lines, which allow us to compare whether or not the black 
hole masses estimated with them are consistent with one an- 
other. As can be noted from column 5 of Table O where BH 
mass estimates based on H/3 lines are in parentheses, the 
differences can be even as large as 0.95 dex (e.g., CSO 879, 
Fig. [2]). Also, we note that in all eight cases the Mg n based 



mass estimate is smaller than the one based on the H/3 line. 
One possible reason could be the uncertainty in the slope, /3 
of the r-L relation; from Eq.[2] /3 = 0.62 ±0.14, so the nom- 
inal uncertainty is ~ 25 per cent. Recent studies on the r-L 
relation also indicate that the slope of the radius-luminosity 
relation is close to f3 — 0.5, consistent with basic predictions 
of photoionization models (e.g., Davidson & Netzer 1979; 
Osterbrock & Ferland 2005). The presence of steeper slopes 
in some analyses appear to be caused by not correcting for 
contributions of the host galaxy to the continuum luminosity 
(Bentz et al. 2006; Bentz et al. 2009). Dietrich et al. (2009) 
have addressed this problem by assuming the slope of the 
r-L relation to be 0.5, as seems to hold for the H/3 line, 
and they then calculated the average scaling constant D; 
following their Eq. (6), viz, 



M bh (HP) = M bh (Mg II) = D 



ALa(3000) 
10 46 ergs -i 



( £W#M(MgIL\2800)\ M 



1000 kms~ 



(3) 



where the M;,f,(H/3) values were computed from their near- 
infrared spectra and the FWHM(MgIIA 2800) values came 
from their optical spectra. They find D — (2.0 ± 0.5) x 
10 8 Mq . Using this empirical constant D in Eq. [3] we recal- 
culated the Mg II -based BH masses; the results are given in 
the sixth column of Table [3] As can be noted from columns 
(5) and (6) and illustrated in Fig. [5] the discrepancy between 
the masses from the two lines for the cases where both can be 
evaluated is reduced, and ranges from 0.05 to 0.36 dex, which 
is within the expected systematic uncertainty of black hole 
masses (around a factor of 4) based on single epoch measure- 
ments (Vestergaard & Peterson 2006). These consistent BH 
mass estimates range from ~ 10 s to > 10 9 M© for our sam- 
ple. We note that the prescriptions we have adopted are not 
unique and that several other prescriptions for black hole 
masses in terms of virial parameters obtained from Ha as 
well as H/3 and Mg II also have been considered. These have 
been recently summarized by McGill et al. (2008) who show 
that the level of agreement we have obtained can be im- 
proved upon somewhat with the consideration of a larger 
number of AGN. 



5 EDDINGTON RATIOS AND BLACK HOLE 
GROWTH TIMES 

We have also estimated the Eddington ratio t = L bo i/L e dd, 
where L bo i is taken as 5.9 x ALa(3000A) and 9.8 x 
Ai A (5100A) for Mg II and H/3 , respectively (McLure & 
Dunlop 2004), and L edd = 1.45 x lO 38 (Mi, h /M )erg s~\ 
assuming a mixture of hydrogen and helium so the mean 
molecular weight is fi — 1.15. The results are given in Fig. [5] 
which shows that distributions of sources showing and not 
showing optical microvariability are not significantly differ- 
ent with respect to either BH mass or I. 

To test the reasonableness of estimated Eddington ra- 
tios, we also computed black hole growth times to compare 
them with the age of the Universe (at the time the QSO is 
observed), by using the following equation (Dietrich et al. 
2009), 
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Figure 6. The upper left and middle panels show the histograms of Eddington ratio L bo i/ L^ dd an d log(Mi,i l /MQ) based on the H/3 line, 
while the corresponding bottom panels show the results based on the Mg II line. The shaded regions correspond to sources with confirmed 
optical microvariability while the non-shaded regions are for those sources for which optical microvariability is not detected. The right 
upper and bottom panels respectively show plots of log(Afi,^ /Mg) versus Lj, D ;. Triangles and circles respectively show sources with and 
without confirmed optical microvariability. 



Table 4. BH growth time estimates (in units of Gyr) from the H/3 line analysis 



object z L bo i/L edd = 1.0 L bo i/L edd obs.^_ age of the 

M(,f,(seed) Mj^seed) Universe 

10 M e 1O 3 M 10 5 M q 10 Mq 10 3 M q 10 s M q [10 9 yr] 



CSO 174 


0.653 


0.83 





63 


0.43 


11.81 


8.94 


6.07 


7.44 


CSO 179 


0.782 


0.82 





62 


0.42 


3.42 


2.58 


1.75 


6.73 


CSO 448 


0.317 


0.74 





54 


0.33 


3.20 


2.33 


1.46 


9.89 


CSO 879 


0.551 


0.79 





59 


0.39 


3.03 


2.26 


1.49 


8.08 


Mrk 1014 


0.163 


0.70 





50 


0.30 


2.49 


1.78 


1.06 


11.43 


PG 0832+251 


0.331 


0.79 





59 


0.39 


6.59 


4.92 


3.25 


9.77 


PG 0923+201 


0.192 


0.84 





64 


0.44 


42.22 


32.19 


22.16 


11.12 


PG 0931+437 


0.458 


0.86 





66 


0.46 


17.29 


13.28 


9.27 


8.73 


PG 1049-005 


0.359 


0.81 





60 


0.40 


5.37 


4.03 


2.69 


9.52 


PG 1259+593 


0.474 


0.77 





57 


0.37 


1.22 


0.90 


0.58 


8.62 


PG 1307+085 


0.154 


0.75 





55 


0.35 


6.83 


5.01 


3.18 


11.53 


PG 1309+355 


0.183 


0.77 





57 


0.37 


5.15 


3.81 


2.47 


11.22 


PG 1444+407 


0.268 


0.75 





55 


0.35 


3.59 


2.63 


1.68 


10.35 


Q 1252+0200 


0.344 


0.93 





73 


0.53 


92.97 


72.91 


52.85 


9.66 


TON 52 


0.434 


0.74 





54 


0.34 


1.61 


1.18 


0.74 


8.91 


US 1867 


0.514 


0.75 





55 


0.35 


1.45 


1.06 


0.68 


8.33 


US 3150 


0.469 


0.76 





56 


0.36 


2.80 


2.06 


1.32 


8.65 


US 3472 


0.532 


0.80 





60 


0.40 


3.98 


2.98 


1.98 


8.20 


US 995 


0.227 


0.75 





55 


0.35 


9.41 


6.91 


4.40 


10.76 



a Values of I estimated using H/3 lines (Table [2}. 
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Figure 7. Observed variation of FWHM (upper) and EW (lower) with Eddington ratio (I = ^bol/^edd) based on both H/3 (left) and 
Mg II (right) lines. Triangles and circles respectively show sources with and without confirmed optical microvariability. The dashed 
lines show the linear regression fits treating £ as the independent variable. The correlations of both FWHM and EW 
with £, suggest that the Eddington ratio is one of the fundamental parameters responsible for some AGN properties. 



M bh (t obs ) = M,T d (io) exp (e ^—^ -?-) , (4) 

\ £ tedd J 

where r = t obs — £rj is the time elapsed since the initial 
time, to, to the observed time, t b s \ M^ ed is the seed BH 
mass; e is the efficiency of converting mass to energy in the 
accretion flow, and t e dd is the Eddington time scale, with 
tedd = a T c/4TrGm p = 3.92 x 10 8 yr (Rees 1984). We used 
Eq. [4] to derive the times, r, necessary to accumulate the 
BH masses listed in Tables [5] and [3] for seed black holes 
with masses of M^ ed = 1OM , 1O 3 M , and 1O 5 M , re- 
spectively. Two cases are considered: (i) BHs are accreting 
at the Eddington-limit, i.e., £ — l^boi/^edd = 1.0 and the effi- 
ciency of converting mass into energy is e = 0.1; (ii) BHs are 
accreting with our observed Eddington ratios and e = 0.1. 
These results are summarized in Tables U and [S] 

In Fig. [7] we show the observed variations of 
FWHM and EW with £ based on both H/3 and 
Mg II lines. There appear to be linear relations of 
both FWHM and EW with £ in these log-log plots. 
For the FWHM plots this is unsurprising since the 
BH masses are proportional to L e dd- To quantify any 
such linear relations for the EW plots we perform 
linear regressions, treating I as the independent vari- 



able, and find 

logEW(H/3) = (1.12 ± 0.14) + (-0.22 ± 0.15)logi, 
logEW(Mg II) = (1.04 ± 0.05) + (-0.34 ± 0.12)logl (5) 

Here the errors on the fit parameters are purely sta- 
tistical. We have calculated the Spearman rank cor- 
relations of logEW with logi?, and found the corre- 
lation coefficient for r 3 (H/3) — —0.22, with null prob- 
ability p nu u — 0.35, and so only a possible weak cor- 
relation is present. Whereas r s (Mg II) = —0.52, with 
Pnuii = 0.01 and so this negative correlation is signif- 
icant. 



6 DISCUSSION AND CONCLUSIONS 

Several recent papers have used the very large numbers of 
quasars discovered by modern surveys to estimate their BH 
masses using virial approaches (e.g., Shen et al. 2008; Fine 
et al. 2008; Vestergaard & Osmer 2009). The large numbers 
of quasars whose spectra are fit in these particular papers 
(between 1100 and almost 57,700) have allowed for very im- 
portant new conclusions to be drawn about quasar demog- 
raphy. Shen et al. (2008) have found that the line widths of 
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Figure 8. As in Fig. [7]with Lbol as the independent variable. This panel indicates that most of the sources with optical microvariability 
are of low luminosities. 



H/3 and Mg n follow lognormal distributions with very weak 
dependencies on redshift and luminosity. Their comparison 
of BH masses of radio-loud quasars (RLQSOs) and BALs 
with those of "ordinary" quasars (RQQSOs) shows that the 
mean of virial masses of RLQSOs is 0.12 dex larger than 
that of ordinary quasars, while that of BALs is indistinguish- 
able from that of ordinary quasars. Fine et al. (2008) used 
Mg II lines to estimate BH masses, and found that the scat- 
ter in measured BH masses is luminosity dependent, showing 
less scatter for more luminous objects. Vestergaard & Osmer 
(2009) have studied the mass functions of BH masses at dif- 
ferent redshifts, and found evidence for cosmic downsizing 
in their cosmic space density distributions. 

The sample we have considered is much smaller, but 
each member has been carefully selected to be among the 
special group of RQQSOs and Seyfert galaxies already ex- 
amined for optical microvariability (e.g., Carini et al. 2007). 
This criterion demands that modest aperture (usually 1— 
2 m) telescopes can make precise photometric measure- 
ments in just a few minutes, and so limits the members 
to the rare QSOs with bright apparent magnitudes (usually 
mv < 17.5). In addition, special care was taken in the fitting 
of the line profiles as discussed in Sections 3.1 and 3.2. 

As the precision of CCD based differential photometry 
has typically improved to better than 0.01 mag for these 
measurements made in a few minutes, the question as to 



the presence of microvariability in RQQSOs has been clearly 
settled in the affirmative (Copal-Krishna et al. 2000, 2003; 
Stalin et al. 2004a,b; Gupta & Joshi 2005; Carini et al. 
2007; Gupta & Yuan 2009; Ramirez et al. 2009). However 
these papers find that the duty cycle (DC), or fraction of 
nights such behavior is detected, is low, roughly 10-20%, 
for RQQSOs observed for reasonably long periods (4 hours 
or more) during a night. On the other hand, blazars show 
much more rapid variability, with DCs for such monitoring 
periods around 50-60% and when observed for > 6 hours 
the blazar DCs rise to 60-85% (e.g., Carini 1990; Sagar et 
al. 2004; Gupta & Joshi 2005). The difference in microvari- 
abilty amplitudes and DCs between RQQSOs and blazars 
can be easily understood if all the fluctuations arise from 
a jet very close to the central engine, but that jet only 
escapes that region and emits significant radio power for 
RL objects (e.g., Copal-Krishna et al. 2003; Gopal-Krishna, 
Mangalam & Wiita 2008). And most RLQSOs that are not 
blazars would only have modest (if any) Doppler boosting, 
so that they also have DCs comparable to those of RQQSOs, 
as found in the samples of Stalin et al. (2005) and Ramirez 
et al. (2009), or at a level between those of RQQSOs and 
blazars (35-40%, in the sample of Gupta & Joshi 2005), 
is not surprising. Only recently has a decent sized sample 
of radio intermediate quasars (RIQs), defined as those rel- 
atively rare objects with radio to optical flux ratios in the 
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object z L bo i/L edd = 1.0 L bol /L edd obs.^ age of the 

M(,h(seed) M(,( l (seed) Universe 







10 M 


1O 3 M 


10 5 M e 


10 M Q 


1O 3 M 


1O 5 M 


[10 9 yr] 


CSO 18 


1.300 


0.83 


0.63 


0.43 


2.59 


1.96 


1.33 


4.73 


CSO 21 


1.190 


0.84 


0.64 


0.44 


3.51 


2.68 


1.84 


5.07 


CSO 233 


2.030 


0.79 


0.59 


0.39 


0.74 


0.55 


0.36 


3.18 


CSO 879 


0.549 


0.75 


0.55 


0.35 


1.34 


0.98 


0.63 


8.09 


PG 0935+416 


1.966 


0.86 


0.65 


0.45 


0.99 


0.76 


0.53 


3.28 


PG 0946+301 


1.220 


0.83 


0.63 


0.43 


1.06 


0.81 


0.55 


4.97 


PG 1206+459 


1.155 


0.85 


0.65 


0.45 


1.01 


0.77 


0.53 


5.18 


PG 1248+401 


1.032 


0.83 


0.63 


0.43 


1.40 


1.06 


0.72 


5.63 


PG 1254+047 


1.018 


0.86 


0.66 


0.46 


3.59 


2.76 


1.92 


5.68 


PG 1259+593 


0.472 


0.80 


0.60 


0.40 


2.68 


2.01 


1.34 


8.63 


PG 1338+416 


1.204 


0.85 


0.65 


0.45 


2.50 


1.91 


1.32 


5.02 


PG 1630+377 


1.478 


0.86 


0.66 


0.46 


1.33 


1.02 


0.71 


4.25 


Q 1628.5+3808 


1.461 


0.81 


0.61 


0.41 


1.36 


1.02 


0.69 


4.30 


Q J0751+2919 


0.912 


0.84 


0.64 


0.44 


1.31 


0.99 


0.68 


6.12 


UM 497 


2.022 


0.82 


0.61 


0.41 


0.78 


0.59 


0.40 


3.19 


US 1420 


1.473 


0.79 


0.59 


0.39 


0.81 


0.61 


0.40 


4.27 


US 1443 


1.564 


0.84 


0.64 


0.44 


1.58 


1.21 


0.83 


4.05 


US 1498 


1.406 


0.80 


0.60 


0.40 


1.05 


0.78 


0.52 


4.44 


US 1867 


0.513 


0.76 


0.56 


0.36 


1.81 


1.33 


0.85 


8.33 


US 3150 


0.467 


0.73 


0.53 


0.33 


1.63 


1.18 


0.74 


8.67 


US 3472 


0.532 


0.79 
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0.38 


2.91 


2.17 


1.42 


8.20 
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0.434 


0.75 


0.55 


0.35 


2.09 


1.53 


0.98 


8.92 


CSO 174 


0.653 


0.82 


0.62 


0.42 


8.21 


6.21 


4.20 


7.43 


TON 34 


1.925 


0.89 


0.69 


0.49 


0.95 


0.73 


0.52 


3.35 


CSO 179 


0.782 


0.82 


0.62 


0.42 


2.92 


2.20 


1.49 


6.73 


PG 1522+101 


1.328 


0.88 


0.68 


0.48 


1.84 


1.42 


1.01 


4.65 



a Values of I estimated using Mg II lines (Table [3)1 . 



regime between those that are truly radio quiet and those 
normally defined as RLQSOs, been targeted for microvari- 
ability monitoring (Goyal et al. 2009); their key result is that 
the RIQs also probably have a DC ^ 20%. Therefore they 
are not likely to be beamed versions of RQQSOs as has been 
frequently suggested (Goyal et al. 2009). 

We conclude that there may be a weak negative cor- 
relation between H/3 EW and the Eddington ratio, I, but 
there is a significant one between the Mg II EW and I 
(Fig. [3 Section 5). This latter point has been made inde- 
pendently and more firmly by Dong et al. (2009b) using a 
larger sample drawn from SDSS. We can also see from Fig. 7 
that there is a decline in FWHM with £; this is unsurprising 
since the BH masses are proportional to L edd . We see from 
Fig. 8 that there is a tendency for sources with detectable 
optical microvariability to have somewhat lower luminosities 
than those with no such detections. This is interesting, and 
should be confirmed by examining larger samples. But it is 
not surprising, as regardless of whether the microvariability 
arises in accretion discs or jets, more massive BHs have nat- 
urally longer characteristic timescales that will presumably 
make any variations harder to detect over the course of a 
single night. 

We also find that the BH masses estimated from the 
FWHMs of both the H/3 and Mg II lines are reasonable, 
in that growth to their estimated masses from even small 
seed BHs is easily possible within the age of the Universe at 
their observed redshift if the mean I values are close to unity 



(Tables 4 and 5). This remains true for the great majority 
of RQQSOs if the value of £ we compute from the current 
continuum flux was constant until the time we observe them; 
however, this assumption does not work for 4 out of the 19 
QSOs with H/3 lines. It is also problematical for 1 out of 
the 26 QSOs with Mg II lines (i.e., CSO 174, which is also 
among the difficult set of 4 using H/3 profiles). 

We find no difference between the EWs (of both the 
H/3 or Mg II lines) and the presence or absence of radio emis- 
sion in those QSOs (Figs. 4 and 7). As discussed in the intro- 
duction, if much of the optical emission in RLQSOs comes 
from a jet, then we would expect the EWs of the RLQSO 
sources to be significantly lower than those of the RQQSOs, 
and if anything, they are slightly higher in our sample. This 
result does not support the hypothesis (e.g., Czerny et al. 
2008) that RQQSOs possess jets that are producing rapid 
variations. Instead it may indicate that variations involving 
the accretion disc (e.g., Wiita 2006) play an important role 
here. 

Improvements to our results could be obtained through 
extensive searches for INOV in a larger sample of RQQSOs. 
It should be useful to divide the sample of sources based on 
their EWs as available from SDSS, or otherwise uniformly 
obtained, spectra. Such samples could be divided into three 
classes based on EW of H/3 , e.g., EW<40A, 40A< EW < 
80A and EW ^ 80A. Such samples should be made as ho- 
mogeneous as possible on the basis of apparent magnitudes, 
z and My. With such larger and homogeneous samples, the 
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absence of a correlation between EW and DC of INOV, as 
found here in our modest sample, could be confirmed or 
shown to be unlikely. 
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